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We derive the Gross-Pitaevskii equation in two-dimension from the first principles of two- 
dimensional scattering theory. Numerical calculation of the condensate wave function shows that 
atoms in a two-dimensional harmonic trap can be condensed into its ground state. Moreover, the 
ground state energy and the wave function of two-dimensional vortex state are also obtained. Quan- 
titative comparisons between 2D results and 3D ones are made in detail. 
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The recent experimental realization of Bose-Einstein 
Condensation (BEC) of dilute alkali Bose atoms in mag- 
netic traps 0^1^] has generated immense interest and ac- 
tivities in theoretical and experimental physics. An inter- 
esting, but much less addressed, question is whether such 
a phase transition due to global coherence exists in two- 
dimensional (2D) space. It has been known that 2D Bose 
atoms do not have a long-range order, and consequently 
BEC of a uniform (untrapped) Bose gas can not occur 
in 2D at a finite temperature since thermal fluctuations 
destabilize the condensate 

There has been no direct experimental observation of 
2D BEC yet, except the recent related works such as 
quasi-condensate of atomic hydrogen |^ and quasi-2D 
gas of laser cooled atoms j^]. Nevertheless, recent theo- 
retical studies have suggested a possibility of BEC in 2D 
in a trapped condition. However, those theoretical ap- 
proaches so far have been only second-handed. Tempere 
and Devreese studied harmonically interacting bosons in 
2D and calculated the critical temperature by ther- 
modynamic arguments to show the possible occurrence 
of 2D BEC. However, it could not provide any informa- 
tion on the 2D-condensate wave function necessary for 
the atomic density profile. 

An analytical approach to interacting bosons in a 2D 
trap by Gross-Pitaevskii equation (CPE), a nonlinear 
Schr'odinger equation for the macroscopic wave function 
of weakly interacting bosons, was also attempted by Gon- 
zalez et. al. but they inevitably faced with a di- 
mensional inconsistency in the relation between coupling 
constant and scattering length. In other words, the well- 
known 3D result, Uq — 47r?i^ a/M (interaction strength 
Uo, s-wave scattering length a, and atomic mass M) is 
not directly applicable to 2D GPE. Bayindir and Tanatar 
recently studied the 2D GPE within the two-fluid model 
of the mean-field many-body quantum statistical theory 
and could calculate the condensate fraction ||l^. How- 
ever, their results were obtained not from the direct so- 



lution of the 2D GPE, but from a density estimation. 
Therefore, the 2D condensate wave function still remains 
unsolved. 

In this Letter, we first derive the "correct" 2D GPE 
from the first principles of two-dimensional scattering 
theory and then show the 2D BEC directly by calcu- 
lating the condensate wave function of trapped atoms. 
The 2D time-independent Schrodinger equation with the 
interaction potential U{p) is given by 



(1) 



where = 2ii„iE/1i'. Here, E is the energy and 
fimi= M/2) is the reduced mass of two identical par- 
ticles. 

The general solution of Eq. ([|) can be obtained in 
terms of the 2D Green's function such that 



Then, the 2D wave function can be expressed as 
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Here Hq is the first-kind Hankel function of order zero 
defined by Ho{x) — Jo{x) -\-iNo{x), where Jq (Nq) is the 
Bessel (Neumann) function of order zero. 

For large p, it shows the asymptotic behavior of 
Hq{x) ~ ^2/7ra; e'(^-"'/4). Note that the Born approxi- 
mation was used and k' = kp. 
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where Fk is the first-order Born scattering-ampHtude in 
2D with a dimension of [length] Note that this asymp- 
totic formula is a solution of Eq. (|l|). Since the scattering 
cross-section in 2D should have the dimension of [length] , 
we can find 
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where ji„ (jsc) is the incident (scattered) flux density. 
Note that the total scattering cross-section for small 
k becomes 27r|Ffc|^ in 2D instead of 47r|/3£)|^ in 3D 
if 3D is the 3D scattering amplitude). Now, if one as- 
sumes a delta- function type interaction such that U{p) — 
{MUo/h?)5'^{p) in Eq. (||), the Fk can be written as the 
following complex form 
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The next step is to find the relation between the scat- 
tering amplitude and scattering length in 2D. Since the 
well-known relation f^D = —a, where a is the s-wave 
scattering length in 3D, is not directly applicable to 2D, 
we need to obtain the relation from a partial wave anal- 
ysis of the two-dimensional collision theory. Outside the 
range of the potential, the scattered wave has a phase 
shift (5m, and the scattering matrix becomes Sm = e^^^"^ 
(note that \Sm\ = 1 for elastic scattering). Therefore, for 
large p, the solution of Eq. (|l|) can be written as 
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where m is an integer. 

The scattering length is usually defined as the distance 
for which the two-body wave function is cut off at zero 
energy. Therefore, the phase shift due to scattering by a 
potential in 2D can be expressed in terms of the scatter- 
ing length so that [O 



1 



2 In kc 



(8) 



where < fca < 1. Note that So ^ ~ka in 3D. Eq. (||) 
can be easily obtained by the following simple argument. 
The scattering length is an equivalent hard-disk radius 
a, imposed by the boundary condition ■0fc(a) = for the 
asymptotic wave- function. For m = 0, the wave function 
can then be expressed as, at large distance and for small 
k 
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Note that a necessary condition for the validity of the 
Born approximation is that the phase shift 5q be very 
small for small k, which can be easily confirmed from 
Eq. (||). Note also that unlike the 3D case where a can 
be negative as well, we do not consider the negative scat- 
tering length since the centrifugal potential of the lowest 
partial wave is negative in 2D [see Eq. (|lj)] so that the 
extrapolated local wave function cuts the radial axis al- 
ways above the origin (i.e., a > 0) [|ll|| . 

Now, from Eqs. (|), (|) and (@), we can express the 
scattering amplitude as a series expansion, and also find 
that m = term is a dominant contribution to the 2D 
system. Therefore, the 2D scattering amplitude becomes 
a complex form, given by 
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Finally, we can obtain the 2D interaction potential C/q 
from Eqs. (|) and (0) as 
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Here, assuming the system is in the ground state of 2D 
harmonic potential, we may approximate the small k to 
be l/auo, where auo — y^h/Muj [uj is the trap frequency). 
Eq. (nil) is a key result of this Letter, and there will 



be no dimensional inconsistency in deriving the 2D GPE. 
Now, the 2D condensate wave function can be obtained 
from the 2D GPE such that 
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where / d'^p4>^{p) — 1 and N is the number of con- 
densate particles. We assume a 2D isotropic, harmonic 
trap potential Vext{p) = ^Mw^p^. Then we can sim- 
plify Eq. ([l^ ) for numerical calculation by introducing 
the dimensionless variables (p — > ahoP, p hujp and 
and the dimensionless interaction of N atoms. 
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where 4>{p) = u{p)/ ^/p. Here the 2D chemical potential p 
can be obtained from the normalization condition. Note 
that the scattering length a is the only atomic parameter 
that contributes to the condensate state. 

The procedure to solve Eq. (|l^ or Eq. ( p^ is 
very similar to the case of 3D |13| . For the non- 
interacting case, the solution is still Gaussian with 4>{p) = 
■K~^l'^e~f 1"^ . In the strongly repulsive limit (Thomas- 
Fermi limit), we obtain the parabolic solution ^^(/c) = 
(2/i — p'^)jU. With a typical value of ka ^ a/ato — 10"'^ 
in Eq. ([ll|), we have plotted the 2D condensate wave 
function versus p and N in FIG. 1. Note that although 
the 2D condensate wave functions show the overall be- 
haviors similar to those of the 3D case, they approach the 
parabolic limit more rapidly with increasing the number 
of atoms. In other words, the effect of atomic interaction 
potential becomes more prominent in 2D. Note also that 
the 2D condensation is not very sensitive to the scattering 
length due to its logarithmic dependence. This implies 
that each bosonic alkali atom with positive scattering 
length may exhibit similar condensate characteristics in 
2D. Refer to TABLE 1 for more detailed comparisons 
between our 2D results and well-known 3D ones. 

The ground-state energy of the 2D system can be ob- 
tained from the energy functional E as n% 
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With the Gaussian trial function <j){x) 
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Eq. (03) satisfies the following rc- 
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The ground-state energy per particle of Eq. ( pii|) is plot- 
ted in FIG. 2 with the same parameter ka ~ a/a/io = 
10~^ and is compared with the well-known 3D result of 
®^Rb. It is interesting to note that as the number of 
atom is increased, the 2D system becomes much unsta- 
ble in contrast to the 3D case. 

Now, let us consider the vortex states in 2D in con- 
nection with superfluidity of the hydrodynamic theory. 
The 2D system may rotate around the 2D trap center, 
resulting in a quantized circulation of atoms. An angular 
momentum quantum number k can be then assigned to 
the quantum winding of the 2D vortex state, which can 
be written as 0(p) = tp{p)e^^'^''\ where V(p) = \/n{p) is 
the modulus. When the phase S is chosen as k9 {k is an 
integer), one finds vortex states with a tangential veloc- 
ity V = nh/Mp. As a result of the quantum circulation. 



the angular momentum of the system with respect to the 
p — axis becomes L — Nnh. Including the vortex term 
in Eq. (|l^), we now obtain the 2D GPE with vortex 



states as 
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The wave functions for k = 1 vortex state are plotted 
in FIG. 3. We observe, in particular, that the overall 
shape and the iV-dependence of the vortex-state wave 
functions are qualitatively similar to those of 3D p2[ . 
Note that 2D BEG transition looks analogous to the 
Kosterlitz-Thouless (KT) vortex-state transition. How- 
ever, the phase transition of 2D BEG does not necessarily 
involve any interaction between atoms and this is the fun- 
damental difference between the two transitions in 2D. 

Although there are fundamental differences between 
our purely 2D results and 3D ones as summarized in TA- 
BLE 1, it will be interesting to consider the quasi-2D 
or (2-|-l)D scheme as a limiting case of 3D, in the sense 
that the atoms obey 2D statistics while their interac- 
tions are 3D. To compare the situation of 2D trap with 
an extremely squeezed 3D trap, we take the 3D external 
potential as Vext{r±,z) = {l/2)Mui^{r']^ + X^z^). Fixing 
the number of atoms and increasing the A, (note that as 
A ^ oo, Vext approaches 2D trap). Solving the 3D GPE 
with fixed number of atoms, we have presented how the 
condensate wave functions behave with A in FIG. 4. We 
can observe the 3D wave functions merge to the 2D limit 
but slowly as A is increased. 

In summary, we have developed the 2D collision the- 
ory to obtain the 2D scattering length and the interaction 
strength. We then derived the 2D nonlinear Schrodinger 
equation (GPE) for trapped neutral Bose atoms with and 
without vortex state. The 2D condensate and vortex 
state wave functions are calculated numerically from the 
2D GPE. Despite the apparent differences, we observe 
the qualitative behaviors similar to the 3D BEG, but do 
not expect 2D BEG for negative scattering length. The 
ground state energy per particle is calculated and com- 
pared with that of the well-known 3D traps. We also have 
considered the practical quasi-2D case as an extreme of 
very thin 3D trap. 
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TABLE I. The comparison of the 3D and 2D BEC. 
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FIG. 3. The 2D wave functions of vortex state for a quantum circulation k=1 and a typical value of feo = 10 ^ . 
0.6 



(})iD 0.3 




FIG. 4. The behavior of asymmetric 3D condensate wave functions for fixed N = 100. From the top 
left A = 10, 100, 300, 1000. The dotted line is the 2D condensate wave function of the same N . Note that 
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